Creating excitonic entanglement in quantum dots through the optical Stark effect 
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We show that two initiaUy non-resonant quantum dots may be brought into resonance by the 
apphcation of a single detuned laser. This allows for control of the inter-dot interactions and the 
generation of highly entangled excitonic states on the picosecond timescale. Along with arbitrary 
single qubit manipulations, this system would be sufficient for the demonstration of a prototype 
excitonic quantum computer. 
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Semiconductor quantum dots (QD's) are often de- 
scribed as 'artificial atoms' due to the discrete energy 
level structure which results from their three-dimensional 
confinement [l|. Consequently, many of the techniques 
of quantum optics are now being used in QD studies and 
have led to the observation of Rabi oscillations ^ , pho- 
ton antibunching Q, and recently the optical Stark ef- 
fect U in single QD's. Such experiments have stimulated 
a great deal of interest in possible applications of QD's in 
quantum information processing (QIP) devices 0,lal3- 

In this paper, we shall analyze the behaviour of two ad- 
jacent self-assembled QD's addressed by an external clas- 
sical laser field, with the aim of controlling the electronic 
interactions between them. We shall demonstrate that it 
is possible to generate and maintain long-lived entangled 
excitonic states in such QD's through the inter-dot res- 
onant (Forster) energy transfer |^ Q . This is achieved 
with a single laser that dynamically Stark shifts the ex- 
citon ground states in and out of resonance, effectively 
switching the inter-dot interaction on and off. 

Our model considers only the ground state (no ex- 
citon) and first excited state (single exciton) in each 
dot, and these two states define our qubit as |0) and 
|1) respectively. Each QD is assumed to be within the 
strong-confinement regime where their typical sizes are 
much smaller than the corresponding bulk exciton ra- 
dius, which is determined by the electron-hole Coulomb 
interaction. As a result, the confinement energy due to 
QD size dominates and mixing of the single-particle elec- 
tron and hole states due to their Coulomb interactions 
may be neglected P]. Any associated energy shift can 
be absorbed into the exciton creation energy; this shift 
is important as it ensures that the resonance condition 
for single-particle tunneling is not the same as that for 
resonant exciton transfer. Additionally, we consider only 
weak inter-dot interaction strengths (~ 0.1 meV) which 
would be expected for two dots with relatively large spac- 
ing (~ 10 nm) [Ifll- Therefore, we may neglect inter-dot 
tunneling of electrons and holes. 

The Hamiltonian for two coupled dots in the presence 



of a single laser of frequency loi may be written in the 
computational basis {|00), |01), |10), |11)} as [Ti = 1): 



H{t) 
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r^i COS ujit 

Q2 COS LUit 

wt + Vxx 

(1) 

where luq is the ground state energy, lui(^2) the exciton 
creation energy for dot 1(2), and lut — luq + lui+lu2. The 
coupling terms Vp and Vxx are the Forster (transition 
dipolc-dipole) and biexciton [y, Q (static dipole-dipole) 
interaction strengths respectively. 

We have assumed that each dot may couple to the 
laser with a different strength, governed by the respec- 
tive Rabi frequency fii or il2, with ili{t) = —2di •E(r, t), 
for i = 1, 2. Here, d^ is the inter-band ground state tran- 
sition dipole moment for dot i, and E(r,t) is the laser 
amplitude at time t and position r. Natural size and 
composition fluctuations in self-assembled dot samples 
(for example in InGaAs QD's [ll|) lead to a large range 
of possible transition dipole moments for each dot. The 
size of the ground-state dipole mismatch (and related ex- 
citon energy difference) between two dots is an important 
factor in determining our ability to control their interac- 
tions. 

We shall first analyze the Hamiltonian of Eq.^ within 
the rotating wave approximation (RWA). This will al- 
low us to derive approximate conditions governing the 
behaviour of the system, and elucidate the mechanism 
for controlling excitonic entanglement. Subsequently, we 
shall characterize the small corrections to these RWA so- 
lutions and perform a full numerical solution of H(t). 

Transforming Eq. ^ into a frame rotating with the 
laser frequency lui with respect to both qubits, we ob- 
tain (within the RWA): 
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where bi = LUi — ui is the detuning of the laser from dot 
i, and an irrelevant constant has been subtracted from 
each term on the diagonal. In order to demonstrate con- 
trol over the interaction V-p we would like to isolate the 
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FIG. 1: Laser induced anticrossing in the {|01), [10)} sub- 
space, for fixed ratio Q.\IQ.-x = 0.55. Eigenvalues are calcu- 
lated from Eq. mwith 5i = 292.59 meV, ^2 = 290.59 meV, 
and Ff = 0.1 meV. 



behaviour of the {|01), |10)} subspace in which it acts. 
We may proceed, utilizing degenerate perturbation the- 
ory, provided that the following conditions are satisfied: 

1^1-^21, |Vf|, |a/2|«mzn(|5,|, |(5, + l/xx|). (3) 

In this case the three subspaces {|01),|10)}, {|00)}, and 
{|11)} are decoupled and we can write an effective Hamil- 
tonian for the {|01), |10)} subspace: 



52 + a^q-li^^ Ff + f^i ^2 (a - /3) 



(4) 



where a = l/(5i, /3 = l/((52 + Vxx), and Q!, = f},/2. 

As we may control both the detunings 8i of the laser 
from the QD's, and the Rabi frequencies fii, it is possi- 
ble to controUably modify the dynamics within this sub- 
space. As is shown in Fig.^ two regimes are of particular 
interest. When the laser is off, the difference in diagonal 
elements can be much larger than the effective interaction 
strength if 



UJX- UJ2~> Vf 



(5) 



In this case, the eigenstates of Eq. 0] approach the com- 
putational basis states |01) and |10), shown away from 
the anticrossing in Fig. ^ that would be expected for 
non-interacting dots. 

In contrast, under the condition: 



br-h=0Jx-uJi^ {^2 - ^f ] (a + /?)> 



(6) 



the diagonal terms of Eq. ^ are equal and the two dots 
Stark shift into resonance under the action of the laser. 
The corresponding eigenstates lie at the anticrossing and 
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FIG. 2: Populations of the four states jOO) , |01) , |11> and |10> 
calculated from Eq. |5| with input |01), pnm = {nm\p\nm), 
Vf =0.1 meV and Vxx = 0. Dashed line: state evolution 
when the laser is always on. The Rabi frequencies and detun- 
ings of the laser are given by f22 = 40.96 meV, fli — 0.55r22 , 
Q,2/'2Si = 0.07, and 5i — 82 = 2 meV. Dotted line: state evo- 
lution when the same laser is on for a time of 7r/(4Kiff) and 
then turned off. The small oscillations in population in the 
7r/(414g) pulse case after the laser is switched off are due to 
some residual coupling between the dots which can be sup- 
pressed by increasing the energy selectivity 5i~S2. Solid line: 
state evolution when the laser is always off. In this case the 
state is almost purely |01) throughout. 



are maximally entangled due to the modified off-diagonal 
interaction 



v,e = VF + n[n'2{c^-P). 



(7) 



They are given by \tJj+) = 2-i/2(|01) + |10)) and |-0_) = 
2-i/2(|io) — |01)). Hence, if the system is initialized in 
the state |01), it will coherently evolve to «|10) during 
the laser pulse, passing through the maximally entangled 
state 2^^/^(|01) + i|10)). This happens with a coherent 
exciton transfer time oi t — tt/{2Vcs). 

Therefore, we may selectively couple the two initially 
non-resonant QD's (satisfying Eq. |S1 before the laser is 
switched on) by the application of a single detuned laser 
satisfying Eq.|Bl which non-adiabatically shifts the eigen- 
states to the anticrossing point in Fig. Q] As soon as we 
wish to decouple the dots again, we simply turn the laser 
off. 

This effect is demonstrated in Fig.|5]where a numerical 
simulation of the evolution of an input state |01) under 
the RWA Hamiltonian of Eq. [21is shown. Without laser 
coupling, the system remains in its initial state |01) with 



a fidelity, T > 1 - 4:{Vf/\Si - 62]) = 0.99. However, 
under the application of a laser satisfying Eq. ^ coher- 
ent oscillations are observed between the states |01) and 
1 10) with an exciton transfer time of 10.9 ps, for the pa- 
rameters chosen here. This transfer is very close to being 
100% complete with little population leaking from the 
{|01),|10)} subspace, justifying our perturbative treat- 
ment. 

This behaviour could be observed in an experiment 
by applying a laser pulse for a series of different pulse 
lengths, and afterwards observing the emitted photons. 
As can be seen in Fig. ^ the two eigenstates when the 
laser is off (which are approximately |10) and |01)) are 
separated by some 2 meV. This is readily resolved in 
a modern spectrometer - and so a measurement of the 
wavelength of the emitted photons can be use to deter- 
mine which of the two dots each one came from. Plot- 
ting the number of each wavelength of detected photons 
as a function of the pulse length would allow a determi- 
nation of the coherent transfer oscillations. Further, a 
t = TT/{4:Vcff) pulse will create and maintain an (approxi- 
mately maximally) entangled state from an initially sep- 
arable one, and this is also shown in Fig|21 Even for the 
small couphng strength (Vp = 0.1 meV) considered here 
this operation is on the picosecond time scale. Therefore, 
we would expect such an entangled state to be long-lived 
in a pair of QD's relative to the timescale of its genera- 
tion. Single and coupled dot exciton lifetimes have been 
measured to be as long as nanoseconds at low tempera- 
tures |12. Ha . [Ijl . Additionally, pure phonon dephasing 
effects are su ppr essed as the temperature is lowered be- 
low 10 K [iiiii. 

We therefore now account for the finite exciton life- 
times by including only spontaneous emission terms in 
the density matrix master equation !la|, and neglecting 
pure dephasing processes: 



<H,p] 



2^ 



Ti (2cr,^ pa+ - cr+cr^ p - pa+a^ 



(8) 
Here, the i label the dipole allowed transitions in the 
coupled system, a^ and a~ are their raising and lower- 
ing operators, and the F^ are their transition rates. These 
terms can lead to a significant reduction of the degree of 
entanglement over time. We can see this by referring to 
Fig. 13 where we show the result of numerical calcula- 
tions which use the full Hamiltonian (Eq. ^ in the mas- 
ter equation (Eq.|HJ|. In Fig.|21(a) the populations of the 
four states are shown as a function of time for the input 
state |01), subject to a square pulse of 5.45 ps duration 
(this satisfies t = tt/{AVcs) for our chosen parameters), 
and subject to significant decay. In Fig. O (b), we plot 
the entanglement of formation (EOF) 17j of the system 
as a function of time for the input state |01), for a va- 
riety of decay rates. The EOF measures the number of 
Bell states required to create the state of interest; for a 
maximally entangled state it is equal to unity while for 
a separable state it is zero. For a general two qubit state 
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FIG. 3: (a) Numerical simulation of the four populations for 
an input state |01) and a 7r/(4Vcff) pulse. The simulation 
uses Eq. |S| which does not rely on the RWA and includes 
exciton decay, but not pure dephasing. The decay rates for 
dots 1 and 2 are given by Fi = t[' — (331 ps)~^ and r2 = 
Tg"^ = (100 ps)~^ respectively. All other parameters are the 
same as for Fig. |5| except Si — 82 = 2.18 meV to account for 
the extra shifts, and lui — 1500 meV. (b) Entanglement of 
formation of the input state |01) as a function of time for a 
7r/(4K:ff) laser pulse, and for a series of different decay rates. 
We keep a constant ratio of ri/r2 ~ |di/d2p — {Q,i/Q,2)'^- 
The calculations are made by using Eas. l8ll51 



it is given by the equation 



Ef{p) = h 



1 + VT" 



(9) 



where h{x) = — a;log2(a;) — (1 — a;) log2(l — x) is the Shan- 
non entropy function, r is the "tangle" or "concurrence" 
squared, which can be computed by using the equation: 



T = C^ = [max{Ai - A2 - A3 - A4, 0}]^ 



(10) 



Here the A's are the square roots of the eigenvalues, in 
decreasing order, of the matrix pp = p a^ 



i<J^P*<J^C. 



where p* denotes the complex conjugation of p in the 



computational basis |00) , |01) , |10) , |11) ^. In Fig. |2| 
we see that at the end of the laser pulse the state has 
become almost maximally entangled and, in the absence 
of decay, stays so once the laser is switched off (the small 
deviation from unity is due to the residual effect of the 
Forster coupling when the laser is off). When decay is 
included the EOF decreases over time, but for typical 
experimental lifetimes of 1 ns |13. Ha . [lj|i it retains a 
value which is higher than 0.94 for times up to 50 ps |19l |. 
The numerical solution of Eq.|Sl required for Fig. 13 has 
been computed without making use of the RWA. The 
behaviour of the system is exactly as expected from the 
RWA case, except for an extra small shift in the dot 
energies when the laser is on. This shift is the largest 
correction to the Stark shift in a perturbative expan- 
sion and arises from the counter-propagating terms in 
the state evolution which are discarded when the RWA 
is made zCr|. The corrections to the RWA can all be 
derived by following the method of Shirley |23 , who con- 
siders the case of a two- level system (which we label and 
1 and consider to have an energy separation of uii) cou- 
pled by an oscillating field (of frequency oji). He shows 
that this kind of problem can be mapped on to a time 
independent one by using Floquet's theory to construct 
a Hamiltonian in an infinite dimensional Hilbert space 
given by 



(mi|HF|/3r7i)=H^^™ 



nujSal3Sn 



(11) 



Here a,/3 G {0, f } and the n,m represent Fourier com- 
ponents of the state evolution. ?i^o is the Fourier com- 
ponent with frequency nuj of the oscillating Hamiltonian 
(which is only non-zero for a ^ (3). Shirley showed that 
the time evolution operator of the two-level system can 
be written as 

Upo.it; to) = J2 (/5"l exp[-i7^F(t - to)] |aO) e™"^*. (f 2) 



He goes on to derive the Bloch-Siegert shift, which is im- 
portant for resonant interactions between the two-level 
system and the oscillating field. However, we are inter- 
ested in the off-resonant behaviour. In the absence of 
any interactions, our system will initially be in a state 
of the form a |00) -I- b |10) (in the notation \an)), and we 
can use the time evolution operator (Eq. I12|l to deter- 
mine the energy separation of the two (eigen)states, |00) 
and 1 10). We follow a similar procedure when the laser 
is switched on. In this case, if the laser is sufficiently de- 
tuned from resonance (i.e. if 5 = (wi —uji) ^ ^/2, with fl 
the laser-system coupling), the \an) are still approximate 
eigenstates and we can employ second order perturbation 
theory to obtain the energy shift in the separation of the 
two levels. We obtain: 



2n 



12 



2n'' 



{S + 2u;i) 



(13) 



where, as before, fi' — ri/2. The extra term represents a 
detuning from resonance by an amount 2uji + S, compared 



to the term with detuning S which is usually the only one 
kept. For the relatively large detunings considered in our 
two dot model, they become non-negligible contributions 
to the Stark shifts for each dot, i, with a magnitude of 
2il[^/{2u!i +6i). Once we account for this extra shift, for 
example in this calculation by redefining the parameters 
Si and ^2 to include it, the system behaves exactly as 
expected from our earlier analysis of Eq. 0] 

We now assess the feasibility of our method by exam- 
ining the state-of-the-art in real systems. In our previous 
work [g, we predicted that the Forster transfer energy 
can be as large as about 1 meV in QD's (corresponding to 
energy transfer times on the sub-picosecond timescale); 
we also showed how a static electric field can be applied 
to the dots to suppress the interaction and so lengthen 
the transfer time as required. In addition to this, exper- 
imental work has suggested that Forster transfer times 
can approach picoseconds in structurally optimized pairs 
of CdSe QD's ,22], and that these transfer times can be on 
the subpicosecond timescale in photosynthetic biomolec- 
ular systems |23 • Molecular systems could therefore pro- 
vide an alternative route to experimental realization of 
the effects we have predicted here. 

QD exciton - laser interaction strengths (which corre- 
spond to the Rabi frequency when the laser is resonant 
with the exciton) of a few meV have been attained in 
an experiment which measured the optical Stark shift 
in AlAs/GaAs hetcrostructures [J. Further, an experi- 
ment which directly observed Rabi oscillations in a In- 
GaAs/GaAs QD photodiode also measured a Rabi fre- 
quency of a similar magnitude |2J|. However, these ex- 
periments were not designed to maximize the laser - exci- 
ton coupling strength, and so the value of 40 meV, which 
we used in our simulations, is not unrealistic. 

The coherent exciton transfer process detailed in Fig. [21 
is equivalent to the realization of iSWAP logi c opera- 
tions between the pair of excitonic qubits |25l |. Along 
with arbitrary single qubit rotations, this gate would be 
sufficient for a demonstration of a prototype excitonic 
quantum computer. Single qubit operations in our two 
qubit system may be achieved by using their frequency 
addressability. A laser resonantly tuned to either of the 
two dots will induce a Rabi oscillation in the dot to which 
it is tuned; in the language of Pauli spin matices this is 
a rotation around the X axis of the Bloch sphere, and 
represents one of the two rotations which arc required 
for arbitrary single qubit operations. The other could be 
obtained in a number of ways. For example a slightly 
detuned laser will cause the qubit state to move around 
another trajectory on the Bloch sphere. Alternatively, a 
higher level \T) (of energy ut relative to the |0)) could 
be used. If this higher level has a different energy for 
each qubit, it could be resonantly excited from, say, the 
|0) state in only one of the two qubits by applying an 
appropriately tuned tt laser pulse. Leaving the system to 
evolve then for a time r before using another tt pulse to 
deexcite it causes the |0) to pick up a phase of ujtt rela- 
tive to the |1). This Z rotation is sufficient to complete 



the set required for universal quantum computing. 

If the ratio of decoherence time to gate operation time 
in our system were large enough, full scale fault toler- 
ant quantum computing (FTQC) J2y| would be possi- 
ble. It is well established that a ratio of around 1000 
is good enough for FTQC, and recent estimates of this 
have been a low as 100 ^3- We have demonstrated that 
an entangling gate can be performed in real systems in 
around 5 ps, which is around 200 times shorter than the 
longest measurements of decoherence times [l^- Hence, 
our proposed scheme could be carried out with a preci- 
sion which is already on the limit required for FTQC; 
incremental improvements in either FTQC protocols or 
in experimental systems should enable the implementa- 
tion of full quantum algorithms. The absolute speed at 
which gates can be carried out in our scheme also makes 
them ideal for smaller scale applications such as quantum 
repeaters i28|, which require much less stringent gate fi- 
delities than full scale FTQC. 

To summarize, we have shown that two non-resonant 
QD's may be brought into resonance by the application 
of a single detuned laser which induces Stark shifts within 



each dot without significant population excitation. This 
in turn allows for control over the inter-dot interactions, 
and hence the generation of highly entangled states on 
the picosecond timescale. The conditions of Eqs. 13 El 
and 1^1 set the upper limits on the energy selectivity and, 
neglecting incoherent processes, the fidelity possible with 
a particular dot sample. The Forster strength Vp sets the 
interaction timescale. In general, as the magnitude and 
difference of the Rabi frequencies increases, and as Vp in- 
creases, so does the feasibility of the proposed idea. We 
believe that the means to demonstrate the effects out- 
lined above already exist. Moreover, these types of ex- 
periments may prove invaluable in assessing the potential 
applicability of semiconductor QD's for future QIP tech- 
nologies. 
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